:arXiv:gr-qc/961 1035v2 22 Nov 1996: 



Quantum state correction of relic gravitons from quantum gravity 



J.L. ROSALES * 

Fakultdt fur Physik, 
Universitat Freiburg, Hermann- Herder-Strasse 3 
D-79104 Freiburg, Germany. 
(February 7, 2008) 



Abstract 



The semiclassical approach to quantum gravity would yield the Schrodinger 
formalism for the wave function of metric perturbations or gravitons 
plus quantum gravity correcting terms in pure gravity; thus, in the 
inflationary scenario, we should expect correcting effects to the relic 
graviton (Zel'dovich) spectrum of the order of (\/m 2 Pl ). These, on 
the other hand, could possibly be measured in a future experiment. 
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I. INTRODUCTION 



In spite of the fact that there still lacks a full quantum gravity theory, a predictive 
semiclassical quantum cosmology has already been developed due to the selection of the initial 
quantum state of the Universe. The latter may be developed upon quantising canonical 
general relativity for the metric tensor corresponding to perturbations about FRW models 
M - [0]. In such an approximation scheme, the parameter representing the relevant energy 
for the initial expansion of the Universe, the effective cosmological constant, A, is much lower 
than the Planck mass Q, 

X/m% < 10- 9 . (1) 

Thus, upon using a Born-Oppenhiemer approximation scheme (that takes into account the 
above mentioned relative small value of matter energy with respect to the Planck mass), 
one can obtain, for a given complex solution of Wheeler-DeWitt equation, the Schrodinger 
formalism for the wave function of metric perturbations |J. Moreover, quantum gravity 
corrections also arise and, at least in principle, they should show how effects of quantum 
gravity might even be measured in the context of these inflationary scenarios. 

The aim of this paper is precisely to show that such effects do exist, in this case, upon 
introducing a shift of quantum gravity origin in the spectrum of the relic gravitons being 
originated during the early inflationary stages of the Universe. Relic gravitons (Zel'dovich) 
spectrum has, as a matter of fact, already been measured from the existing 3K cosmic back- 
ground radiation anisotropy 0; on the other hand, a certain quantum gravity modification 
could be possibly measured in the next future by some experimental devices ||. 

II. SEMICLASSICAL WAVE FUNCTION IN DE SITTER SCENARIO 

The simplest (although still realistic) model for the early stage of the Universe consists of 
a massive scalar field with matter potential U(<p) in a FRW spacetime described by a single 
degree of freedom, the scale factor, a(t). After the Hartle-Hawking prescription [K|, the wave 
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function will depend on the value [/(0(E)) ( E being the boundary of the compact euclidean 
4 - sphere) of the inflaton (scalar field) potential. This is just de Sitter case with U(4>) ~ A 
playing the role of an effective cosmological constant for large initial values of 0; thus, in 
this model, we only take into account, for the wave function, an implicit dependence on the 
inflaton field (in classical terms: we will ignore the backreaction of matter corresponding to 
the kinetic energy of the field, so that matter couples directly to curvature using the effective 
cosmological constant). 

Moreover, the general approach is only physically consistent if we obtain a description 
for the Schrodinger evolution of metric perturbations or gravitons, i.e., the tensor harmonics 
of the three sphere d n . Now, upon defining a = log(a) as the gravitational variable, the 
metric perturbation Hamiltonian is 

^ = E^ + (- 2 -i)e 4Q C (2) 

where, as a consequence of the perturbative character of the {d n }, 

d 2 J\ « 1 (3) 

Therefore, upon assuming U(4>) ~ A, and neglecting second derivatives with respect to the 
matter field, Wheeler-DeWitt equation reads @ 

' 1 +^(-e 4 « + Ae 6a )+tf m ]* = (4) 



l 2m 2 Pl da 2 2 

In the semiclassical approach we write the wave functional *&[a, {d n }} as 

* = e iS (5) 

then we expand 5* in the form 

S = m 2 Pl S + St + mp 2 S 2 + .... (6) 

We now insert the expansion defined by Eqs. (5) and (6) in Eq.(4) and compare expres- 
sions with the same order in mp\. The highest order yields to 



E(^) 2 = o (7) 

n ULl n 

Thus, So depends only on a (the three-metric). 

The next order leads to the Hamilton- Jacobi equation for gravity alone: 

_ (^°)2 + xe 6a - e 4Q = 0. (8) 
da 

If we define a functional ip according to (Si = J2n Sin) 

tfj = D(a)e i ^ s ^ = U n ^j 0n (9) 

where, 

D(af = (^) = S> (10) 

then order mP Pl will imply that ipo is the solution of 

.8 S dijjo _ .di> 



— i- 



= i-£r = H m ifj - (11) 



da da dt 

Which is the functional Schrodinger equation for matter fields propagating on a fixed curved 
background; t is usually called WKB time. 

If we do not take into account the value of S2 in the formal expansion, we get to this 
order of approximation the WKB wave function given by 

^r(l) = _ e im^5 ^ o _ ( 12 ) 

III. QUANTUM GRAVITY CORRECTIONS 

Upon defining S2 = a(a) + r]({d n },a), we obtain the second order correction for the 
wave function 

where a and 77, representing quantum gravity correcting terms to the WKB wave function, 
satisfy, after some algebra f| 



Q/ , _ i ^ d 2 r] t i ^ drj dipp t 1 , „ 1 , „ 

~ 2 ^ ddl %oV 94 94 + 2S^ %b ° 2^ % 1 j 

a= -4^ + 8^ (15) 

Eq. (14) involves entanglement of the modes and, therefore, it is difficult to solve. 
Moreover, in general, field modes and gravitational degrees of freedom are also entangled 
so we do not expect to obtain its general solution. Nonetheless, let us define r]({d n },a) = 
J2nVn(d n ,a). In that case, we obtain, 

^ " ddl dd n dd n + 2S^0„ °" ° 2^0n^°" 1 + 

+ ^E- J — ^'oi (16) 



2 g ^< 



0/ 



where ipon is the solution of the Schrodinger equation. The non boundary proposal B 
picks up the ground state (these results are desired in the semiclassical approach to gravity 

!), 

iP 0n (a,d n ) = N n (a)e- 1 * na2d * (17) 

and as a result of this, last term in Eq. (16) is 0{a i d\d 2 ) which, after Eq. (3) should be 
neglegible for a 2 — > 1/A. Thus the modes are exactly disentangled after the selection of 
this particular initial condition and we can finally write, for the correcting phase of a single 
mode vave function, 

_ ql / _ i^Jh i dr] n dip 0n 1 ,, q/l I ,„ 

boVn ~2dd 2 n + ^ n dd n dd n + 2^o/°" ° 2^ 0r > 

We must now solve the system of Eqs. (15) and (18) in order to obtain predictions from 
the corrected wave fuction. To this aim, it is better to consider the very early stages of 
the Universe where such terms should be relevant, i.e., we must restrict our calculation to 
the region where A is really a constant; now, for a 2 — > 1/A, we can opperate the expressions 
a little bit further. First, since the theory should not depend on the selection of the time 



parameter, we are allowed to make our predictions for the particular semiclassical evolution 
parameter given by 

9(a) = XD(a) = Xa 2 (Xa 2 - 1) 1/2 (19) 
thus, for a 2 = 1/A(1 + 9 2 + 0(9 4 )), we get, after Eqs. (18) and (15) 



i drj n 1 d 2 r] n 1 drj n dip 0n i d 2 ip 0n 2 dip t 



1 nan./. L /i o/i J V^"/ 



A 90 2 ip 0n dd n dd n 29 2 ip 0n d9 2 9 89 

and, 

da _ 5A 
89 ~W 



(21) 



Which shows that there seems to exits an apparent divergency for \l/( 2 ) as 9 — ► 0. Moreover, 
from Eqs. (9) and (17), upon factorizing Van Vleck determinant, we pick up the n - mode 
wave function ip 0n as 

If we now separate the factor ordering dependent part (i.e., that arising from the Van Vleck 
determinant) of t] n in the form 

Vl (9, dj) = Vji(0, dj) + V2(9)5 nj (23) 

we finally obtain, replacing Eqs. (22) and (23) in Eq. (20) 

^[a(9) + V2 (9)] = (24) 

i drj nl ld 2 7] nl 2 {a \ dr lni 

A^ = 2^r-^ na( ^ (25) 

Eq. (24) demostrates that, after the adiabatic approximation, the divergencies arising in a 
and rj n as we approach 9 —>■ cancel out each other exactly. This, on the other hand, is a 
consequence of the fact that the phase correcting terms from quantum gravity should not 
depend on the selection of the factor ordering ||. Therefore, in \E f ( 2 - ) , the only physically 
relevant quantum gravity correcting phase factor is rj n \{9, d n ). 



The initial state of the Universe is taken on the three-sphere a 2 — > 1/A, in this case, we 
can obtain an exact solution of Eq. (25) upon making the obvious substitutions 

Vm = e l£Xe y(d n ), (26) 

for some unknown constant e. Then we write 

d n = (\/ny/ 2 x n , (27) 

k = eX/n, (28) 



using these expresions, Eq. (25) transforms into 

^| - 2x n -p- + 2ky = 0, (29) 

which is the Hermite equation. The requirement of normalizability in d n of the corrected 
wave function ip n = ipo n e m p' imposes that 77^1 should only be given in terms of polynomials 
in d n , i.e., we select the constat k in Eq. (29) 

A; = 0,1, 2, (30) 



or 



Vm(0,d n ) = { \^e^ 2 [2qy/ 2 d n X}+g 2 e^[(^ - 2)A]. (31) 

After Eq. (22), we notice that gi 7^ in Eq. (31) (i.e., a linear term in the phase of the 
n-mode wave fuction of gravitons) means that the expected initial number of gravitons is 
different from zero and it is given by 

N= Jim 1 2 ^ ( : /A D 2 =I#I 2 - < 32 > 
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IV. THE SPECTRUM OF RELIC GRAVITONS 



Since the initial number of gravitons is different from zero, it will also change the mea- 
surable properties of the resulting spectrum corresponding to the statistic of gravitons in 
inflationary models; here, in order to calculate the expected changes, we follow the method 



In the first approximation, the graviton spectrum produced by an inflationary stage is 
enterely independent of the mechanism that produces the inflation. The only inputs which 
are needed to find the graviton spectrum are the classical metric of space-time, and the 
initial quantum state of the gravitational perturbations. 

For convenience, one may assume that the Universe is approximately flat, so that the 
metric takes the form 



In any case, since we will only study the situation corresponding to wavelengths which are 
shorter than the present-day horizon scale, it could be taken as a good approximation if the 
Universe were either spatially open or closed. 

The classical spacetime begins as de Sitter space but then undergoes an instantaneous 
phase transition at, say, t = t\, after which it evolves as a radiation-dominated model until 
the time t = to- The scale factor describing this model is 



of B. Allen jTT . 




(33) 



a(t) = {t/ti)a{ti) for t x < t < t 




(34) 



We would also require the solution of Einstein equations, imposing 



8irGp 
3 




1 



A 



(35) 



Let us now determine the gravitational-particle production in this spacetime. 



A gravitational perturbation with comoving wave number k is represented by 
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K v = a(t) 2 e^(k)(f)(t)e ikx + cc (36) 

where is the polarization tensor. The physical frequency of the wave is u = k/a(t). The 
amplitude (f> also obeys the perturbed Einstein equations, it leads to 

<j> + (2d/a)0 + k 2 (p = (37) 

The choice of a solution to the equation for corresponds to the choice of an initial quantum 
state for the gravitational field. In de Sitter stage, the solution representing a de Sitter- 
invariant vacuum state is 

Mt) = [a(*0/a(t)]{l + t(Xf 2 (a(t)/k)}e~ tk ^ (38) 

On the other hand, if, after Eq. (32), the initial state is given by a small iV-graviton coherent 
state, we must correct the vacuum state upon adding a negative frequency solution to the 
above expression, 

<f>(t) = <f> v (t) + (igi\/m 2 Pl )cf>* v (t) + 0(X/m 2 Pl ) 2 (39) 

The corresponding solution to the wave equation in the radiation stage is 

<f> r (t) = (a(*i)/a(t)){a r e- <fc (*-* 1 ) + (3 r e +ik ^} (40) 

where a r and j3 r are Bogolubov coefficients. 

By matching the modes at t — t\ one obtains for the negative frequency coefficient, 
representing particle creation 

= (-i 9l X/m 2 Pl )(l + ~ zkHj (41) 

On the other hand, upon taking into account that k = a(to)u, t\ = l/a(ti)A 1 / 2 we finally 
obtain, for the total number of gravitons at time to, 

4 a(to) uj* 

where, 
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, 1 - 2(A/<)[2i fe[ffl ] [ |il]^ + / mW (l - 2[4g]"$ (43) 

Therefore, the spectrum should be given by the corresponding energy density dp g = P{u)du 
for a density of states dN = uj 2 duj / (2ir 2 ), i.e., 

dp g = P(u)du = 2tu^-^\p r \ 2 (44) 

or, 

Here, e(u) differs from Zel'dovich's spectrum due to the presence of quantum gravity effects. 



If we take into account some phenomenological values of ~ 10 28 (see also |IH), we 

conclude that such effects would already be present for those frequencies of the order of 
uj ~ 10 13 Hz. 

The result in Eq. (45) may also be expressed in terms of the effective Hubble parameter 
during inflation, H = A 1//2 , and the recombination density pn = Po[^y] 4 , then, Eq. (35) 
taken into account, we obtain the typical value 

UJ dp„ 2 . H ,n / \ / 

Qgraviton = "J^ = — ( ) e(w) (46) 

Pr dw Sir mpi 

where e(u) ~ 1 + 0(\/m Pl ), represents a perturbation to the predicted Zel'dovich plateau. 



V. CONCLUDING REMARKS 

Cosmic Background Radiation is used in order to test phenomenological models of the 
early Universe. In these scenarios, we have just shown that, quantum gravity corrections 
for the spectrum of gravitons, would possibly lie on the range of frequencies technically 
accessible in some future experimental devices measuring CBR anisotropy ||. 

Semiclassical gravity, represented by ty( 2 \ is, in this framework, a predictive and testable 
theory of initial conditions. On the other hand, we have also seen that the adiabatic approx- 
imation for the wave function of metric perturbations, leading to the ground state for the 
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wave function of gravitons, can be thought as being technically correct since, in this case, 
there would not exist formal divergencies depending on the factor ordering for the resulting 
quantum gravity corrections on the initial three-sphere. 
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